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USE OF NONLINEARITIES TO CCMPENSATE FOR TBE EFBECTS OF 
A RATE— LIMITED SERVO ON THE RESPONSE OF AN 
AUTCMATICALIY CONmOLLED AIRCRAFT 
By Stanley F. Schmidt and Nilllam C. Triplett 

SUMMARY 


Hydraulic servos of the type normally used in airplane-autopilot 
ccmhlnatlons are nonlinear over the greater portion of their operating 
ranges. One of the most important nonlinearities is the limit on output 
rate which often results in oscillatory or unstable airplane responses 
to large input commands and sluggish responses to small error signals. 
These undesirable effects can generally he compensated for hy introducing 
other nonlinear elements into the system. 

This report describes a simple method for determining the nonlinear 
gains reqxilred to give optimm responses for step inputs of all magni- 
tudes, This method is based on the fact that the control surface moves 
at its maximum rate during practically the entire transient maneuver and 
thus the servo system can be considered as a simple "on-off" type con- 
troller. The method requires a knowledge of the transfer function that 
describes the airplane response but ignores the dynamics of the servo 
and requires only simple hand calculations. 


INTRODUCTION 


The present trend toward the elimination of the human pilot as the 
primary controller in modern high-speed aircraft has led to severe require- 
ments for autopilot performance. In many applications the autopilot miist 
be capable of controlling the aircraft through violent maneuvers and at 
the same time must give rapid and precise response to small error signals. 

Many of the problems encountered in automatic control of aircraft 
do not appear in other applications of servomechanisms. For example j the 
airplane has a complicated dynamic response that varies over its operat- 
ing range. Furthermore, the autopilot is required to develop large 
forces or moments and still meet strict requirements with regard to size 
and weight. For this reason hydraulic servos are generally employed, 
and oil flow to the cylinder is restricted not only by the limits on size 
and weight but also by structural considerations. The net result, there- 
fore, is a servo with limited output rate. 
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Thus, even though its ave]>-all response characteristics are extremely 
complex, the autopilot for many purposes can he considered as a simple 
linear system with a limit on the rate of control surface deflection. It 
is this limit on the control— surface rate that is often a major source of 
difficulty in designing a system that is stable for large inputs and suf- 
ficiently responsive to very small commands- Adjusting the system parame- 
ters to give a stable response for large Inputs usually results in slug- 
gish responses to small Inputs. 

One way to compensate for such undesirable effects is to design 
system parameters as functions of error. Considerable effort has been 
expended on this general problem of developing nonlinear elements to 
improve the performance of linear as well as nonlinear servomechanisms. 
References 1 and 2 are typical examples in which phase plane methods of 
analysis are used. In this report, however, a different approach was 
necessary in considering the more complex case that is characteristic of 
an airplane— autopilot combination. A method is developed f or ' des ignlng 
appropriate nonlinear functions of error into a rate— limited system to 
give large gain levels for small errors and low gains for large errors so 
that satisfactory responses may be obtained with step inputs of any mag- 
nitude. The method is illustrated by two examples which consider hypo- 
thetical airplane— autopilot combinations with control— surface rate 
limiting. 


NOTATION 


9 pitch angle, deg 

&± pitch— angle input, deg 

cp roll angle, deg 

cpi roll-angle input, deg 

e error, deg 

p Laplace operator 

t time^ sec 

tg control— surface switching time, sec 

Vg servo error, volts 

aileron angle, deg 
6 g elevator angle, deg 

gain constants, volts/deg 
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Kg,K^ gain constants, volts /deg/sec 
Ke(e) nonlinear gain function of e, volts/deg 

K^e) nonlineeir gain function of e, volts/deg/sec 

A dot over symbol Indicates derivative with, respect to time, 

. EFFECTS OF SERVO RATE LIMITIIIG ON TYPICAL 
AIRPLANE AUTOPILOT SYSTEM 


In figiire 1 is shown a block diagram of the roll-control channel of 
a typical airplane— autopilot combination. The system consists of a roll- 
rate feedback for stabilization and a roll— angle signal as an attitude 
reference- The servo is assvmed to have the open-loop transfer function 

^ ^.50 

ve P 


With aileron position feedback, the closed— loop servo response is 


1 

1 + 0.02p 


This simplified representation of the servo is usually justified because 
its response is so much faster than that of the airplane. For this 
example the airplane roll— angle response is defined as 


'P = 8.1 

8a p(l + 0.3p) 


Cross— coupling terms have been ignored and the airplane is assumed to 
have a single degree of freedom in roll. The constants were chosen as 
typical of a high-speed airplane. 

By means of an electronic simulator the gain constants and 

were adjusted to give the' desired response to a step input (a rea- 
sonably fast response with little overshoot) such as shown by the solid 
line of figure 2. Values of the gain constants for this condition are 
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Ke = 3.33 


and 


= 0.41? 


Also shown for comparison is the more sluggish response obtained when 
Kg is reduced to I.67. 

When an aileron rate limit of 50° per second is introduced, the 
solid line of figure 2 is, strictly speaking, descriptive only of respon- 
ses to Inputs of no greater than 0.3° since the maximum rate is attained 
when Ve is 1 volt (fig. 1) . For larger commands the system is non- 
linear as indicated in figure 3. Here aileron- Mod roll— angle responses 
for step inputs of 2°, 5°^ 7*5°> 10°, and 15° are plotted. For an input 
of 2° the aileron rate is limited for only a short time, and the response 
differs very little from that shown in figure 2. For larger inputs, how- 
ever, the response becomes oscillatory and finally unstable for a step 
command of 15°; in eanh case the control surface moves at its maximum 
rate (±50° per sec) until the roll angle has essentially stabilized at 
the command input. By reducing the gain Kg the stability for large 
inputs could be improved, but then the respoixse to small inputs would be 
slviggish, as indicated in figme 2. 

To generalize on the control motions shown in figure 3 > it can be 
stated that to correct an initial error in the shortest time, the control 
surface should move at its maximum rate, reverse direction at precisely 
the right time, and travel at its maximum rate in the opposite direction. 
The fact that a response is sluggish indicates that the control surface 
has changed direction too soon; an oscillatory response means that the . 
surface has traveled too far before changing direction. Thus, with a 
fixed control— s-urf ace rate there is an optimum reversal point, or "switch- 
ing time" tg for each input magnitude. 

During practically the entire transient response (<P;j^> 2.0°), the 
operation of the sjrstem is essentially the same as an "on-off" or "bang- 
bang" type of controller commanding a fixed plus— or-nninus control— s'urf ace 
rate. The only variable quantity is the time at which the rate changes 
sign. 


It is apparent from figure 3 that the proper control— surface motion 
can be attained for step inputs of different magnitudes only if the sig- 
nal to the servo is modified in some fashion. The method described in 
the following section makes use of this fact in establishing gain level 
as a suitable nonlinear function of error. 
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DESCRIPTION OF METHOD 


The proposed technltiue for determining the proper nonlinear gains 
for compensating the effects of rate limiting can "best he described hy 
use of ill\xstratlve examples. The first example considers the system of 
figure 1, while the second example shows how the same basic method can 
be applied to a more complicated situation. 


Example 1 


The first step is to determine the optimum switching time, tg, as 
a function of cp^. This is done by calculating the airplane roll respon- 
ses to a n\miber of constant— rate (50° per sec) triangular Inputs as shown 
in figure 1)-. These may be calculated readily by evaluating the response 
to a constant— rate aileron input and using the principle of superposition. 
Curve (2), for example, shows that If the aileron moves at its maximum 
rate for 0.35 second and then changes direction, the airplane will attain 
a maximum roll angle of 39-^° approximately 0.92 second. Thus, to 
obtain the quickest response to a step input command of 39.^° without 
overshoot, the system gain should be adjusted so that the aileron will 
reverse direction at 0.35 second. 

The values of tg corresponding to the various peak values of cp 
from figure 4 are plotted as the solid line in figure 5- This curve is 
labeled "optimum" because it defines the response requiring the least 
time to reach and remain at zero error. Furthermore, it marks the bound- 
ary between an oscillatory and a no-overshoot response. For purposes of 
comparison, switching times indicated by figure 3 (for Kg of 3-33) ane 
also shown. As expected, the switching times are greater than optimum 
over most of the range, indicating an oscillatory system. When is 

reduced to 1.67^ figure 5 indicates the system to be slTigglsh for Inputs 
less than 20° but still oscillatory for larger commands. 

After determining the desired switching time as a function of input 
magnitude, the next step is to find a nonlinear function that will cause 
the aileron to reverse direction at the proper time. While there are 
several possible choices, the most obvious is the replacement of Kg 
with a nonlinear function of error. At the instant the aileron reverses 
its direction of motion the signal Vg to the servo has just reached 
zero. Thus at time tg (from fig. l) with = 0.4l7 

Ve = Kg€ — 0.4179 — 8 a = 0 (l) 

For a given cp^^, t is taken from figure 5j 9 is calculated as shown 
in figure 6j and Oa is equal to 50 . As outlined in table I the 

desired value of KgS for each <pi is obtained directly from equation 
( 1 ). The error angle e at tg is simply 91 — 9 . 
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In figure "J, Kg, designated as Kg(e), is plotted as a function of 
error and is the gain level necessary for an optimum response. At zero 
error the curve indicates an infinite gain; however, this value is phys- 
ically impractical and has no significance because for small errors the 
system operates primarily in its linear range. Thus with K<p fixed at 
0.4l7» Kg should he restricted to a value of 3-33 shown by the dot^ 
ted line so that in its linear range the system will have the response 
characterized by figure 2. 

By use of an electronic analog computer the response of the system 
with this nonlinear gain function was detenained and is shown in figure 
8. It can be seen that the system has a rapid and stable response for 
inputs as large as 60°. It should also be noted that for the larger 
Inputs, two reversal points occur before the system remains within its 
linear range. While only the first point was considered in the calcur- 
latlon of the nonlinear function, the analog computer results show that 
the succeeding reversal points (which occur whenever equation (l) is 
satisfied) are properly timed for a near optimum response. 

The same improvement in system response may be obtained by making 
the damping parameter K,p rather than Kg the nonlinear function of 
error. The calculations are similar to those for the previous case and 
are also shown in table I. In this case Kg has a fixed value of 3»33# 
and at time tg 


3.336 -K^«|,- 5^ = 0 


The term K^ is plotted as a nonlinear function of e in figure 9 which 
indicates an infinitely large negative gain at zero error. As in figure 
7 this value is not significant and for small errors Km(e) may have the 
value of 0.i|-17 shown by the dotted line. The corresponding system respon- 
ses shown in figure 10 are almost identical to those shown in figure 8, 
Thus, either Kg or Km may be replaced by an appropriate nonlinear 
function of error to obtain satisfactory airplane responses to a wide 
range of step command inputs . 


Example 2 


This example shows how the basic method may be applied to the more 
complicated system shown in figure 11 which is the block diagram of the 
pitch channel of a typical airplane— autopilot combination. 

The transfer function of the airplane contains an exceptionally 
large lead term in the numerator (1 + 2.2p), and hence it was necessary 
to include -the compensating network 
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1 + 


1.2 

1 + 2.2p 


in tlie system in order to oltain a satisfactory response even for tlie 
linear case with no rate limiting. The system parameters were adjiisted 
to give the response shown in figure 12. In this case 


Kg = 10 

H = 1-3 


With an elevator rate limit of 50° per second the response became 
unstable for input commands greater than 2.5° as shown in figure 13 . 

The calculations involved in expressing Kg as an optimum nonlin- 
ear function of error are the same as described in Example 1, except 
that at time tg, when the signal to the servo is zero, the following 
condition applies 

ve = Kg€ (l + 1 - Se = 0 


or 


KgS + 


1.2Kge 

1 + 2.2p 


1.30 + 8e 


( 2 ) 


For each given 0^, the quantities 9, 9, e, and 5g are calculated at 
the proper tg as Illustrated in Example 1 (figs, k to 6 and table l)'. 
Equation (2) may then be written as 


Kg€ + Vj, — 1.30 + Sg 


(3) 


where 


vx 


1.2Kgg 

1 + 2.2p 


(4) 


To solve equation (2) for KgS it is first necessary to evaluate vi at 
the time tg corresponding to each given value of 0^^ . This may be done 
by writing equation (4) in the fonn. 


* A " -iti 
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The solution of this differential equation (as shown for example in ref. 
3 ) with zero initial conditions is 


Vl(tj 


12 

1 ) = e 2.2^ e2.2 Kee(t)dt 


(5) 


However, at this point in the calculations Kg is an tinknown function 
of error and an exact solution is impossible. To obtain a first approx- 
imation of Vi(tg) it may he assumed that Kg remains constant until 
time tg . Thus 


_ tg t 

Vi(tg) ft Kge 2.2 J e2.2 e(t)dt (6) 


Here e(t) = d± — 9{t) where 0(t) is the known response to a constant- 
rate elevator input. This integral can be evaluated conveniently by 
means of graphical or numerical integration to give Vi(tg) expressed in 
terms of Kg. This approximate value- is then substituted into equation 
( 3 ) to give a close if not exact value of Kg for each 9j_. 

More precise results could be obtained by repeating the solution of 
equation (5) with actual values of Kg obtained in the first trial. 

While there may be cases where more than one iteration is required, it 
was found in the present example that vx was small compared to KgS 
and that a second solution was not necessary. The quantity Kg as a 
nonlinear function of error is plotted in figure l4 and the corresponding 
system responses eire shown in figure 15 . Comparison with figure 13 shows 
a marked improvement in the performance of the system. 


DISCUSSION 


The method illustrated in this report was designed to give an opti- 
mum response with zero overshoot. By neglecting the lag of the servo it 
was assumed that the control surface responded -instantaneously at a fixed 
rate when the servo error signal changed sign. "However, when the systems 
were simulated on the analog computer, a representative value of servo 
lag was included. The small overshoots apparent in figures 8, 10, and 
15 are the results of this lag and also of possible inaccuracies in the 
analog simulation. In any event, the results Indicate that the simpli- 
fied method is valid when considering servos with reasonably small time 
constants . 
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In some cases It may te desirable to piirposely design a system to 
have a given overshoot for all input magnitudes. The same basic method 
is applicable; and in Example 1, if a 10-percent overshoot is desired, 
the abscissa scale of figure 5 is merely multiplied by O. 90 . With refer- 
ence to curve of figure k, a switching time of 0.35 second would now 
correspond to an input command of 0.9° of 39-^° or 35-5°» but the response 
would actually reach a peak value of 39 

It is also possible to expand the basic method to Include position 
limits in addition to rate limits. The procedure is the same except that 
the triangular inputs, such as shown in figure 4, are out off at the 
value corresponding to the position limit. The switching time is still 
the point at which the control— surface rate becomes negative. For the 
two examples in this report it was fotmd that the addition of position 
limits would have a stabilizing influence. 

In general, the effects of limiting on a particular system can be 
shown clearly by plots similar to figure 5- These can be used to good 
advantage to obtain qualitative and even rough quantitative measures of 
system performance without resorting to the calculation of ccanplete time 
responses. After determining the optimum curve for a particular system, 
it is a simple matter to plot corresponding ctirves for various fixed 
values of system gain. As long as the switching time remains below the 
optimum curve, the step response has no overshoot (if none exists in the 
linear system with no limiting) . If tg is greater than optimum for a 
particular input, the response is oscillatory or even unstable, necessitat- 
ing a decrease in system gain (Kg) or an increase in feedback gain (K<p) . 

In this regard it is Interesting to note that figure 9 indicates 
that with the nonlinear (rate— limited) system the damping parameter (K^) 
must increase with Increasing error in order to obtain an optimum 
response. This is in direct contrast to a common practice for improving 
the response of linear second— order systems (ref. l) where it is neces- 
sary for the damping to decrease with Increasing error, thus allowing 
fast response to large errors while effectively preventing overshoot. 

In Example 1, the rolling acceleration 9 is proportional in the 
steady state to the control— surface rate 6^. Thus, a limit on Sg, 
effectively limits the second derivative of the output 9, and may be 
termed an output ^acceleration limit. Similarly, a 5g position limit 
would restrict 9 and could be specified as an output-velocity limit. 

In this sense an acceleration limit is generally destabilizing, but for 
the type of systems considered the addition of a velocity limit tends to 
improve the stability. 



10 


MCA TN 3387 


As an extension to the present study a more complete generalization 
of the effects of acceleration and velocity limits would he of Interest, 
This broader investigation should consider systems of varying degrees of 
complexity and also the effects of external disturbances other than pure 
step commands. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Oct. 15, 195^. 
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Servo Airplane 



Figiire 1 .— Block dlagrean of roll-control syBtem used, in Example 1 . 
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Figure U.— Roll— angle responses to constant— rate ( 50 ° per sec) 

aileron inputs . 
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Figure 6.— Rolling— velocity response to a constant rate (50° per sec) 

aileron input . 
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Figure 8.— Response of modified roll— control system to step commands with 
Kg as a nonline6LT function of. error. 
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Figure 10.— Response of modified roll-control system to step commands 
with K(p as a nonlinear function of error. 
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Figure 13 .— Response of 1381810 pltchr-control system to step commands with 

elevator rate limit of 50 ° psr sec. 
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